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Abstract 


The  motion  of  a  freely  floating  or  submerged  body,  which 
is  moving  with  a  constant  average  forward  speed  and  oscillating 
arbitrarily  in  any  of  the  six  degrees  of  freedom,  is  formulated 
as  an  initial-value  problem.  The  seaway  is  assumed  to  be 
arbitrary.  The  body  is  assumed  to  be  'smooth',  but  no  symmetry 
of  the  body  is  required.  The  fundamental  assumption  is  that 
both  the  free- surface  disturbance  due  to  forward  motion  of  the 
body  and  the  oscillations  are  small  enough  so  that  the  problem 
may  be  linearized.  By  an  approach  similar  to  that  of  Wehausen 
(1965),  it  is  shown  how  the  present  treatment  of  the  problem 
leads  also  to  Ogilvie's  (1965)  modified  results  of  Cummins' 
(1962)  decomposition  of  the  velocity  potential  for  the  case  of 
an  oscillating  body  with  a  constant  average  forward  speed.  The 
linearized  equations  of  motion  of  the  body  are  then  derived  as 
a  set  of  six  lntegro-differential  equations.  Existence  and 
uniqueness  theorems  are  not  established  either  for  the  boundary- 
value  problem  or  for  the  integral  equation  which  is  constructed. 


i 


r ' 

i 

f 


i 


.jt 


_  .j 


I .  Introduction 


In  the  study  of  the  problem  of  ship  motions,  It  Is  desirable 
to  be  able  to  write  down  equations  of  motion  which  are  valid 
whatever  the  nature  of  the  seaway.  This  means  that  the  validity 
of  these  equations  should  not  require  the  forcing  functions  to 
depend  sinusoidally  upon  time.  Cummins  (1962)  made  an  important 
advance  toward  this  goal  by  considering  a  certain  decomposition 
of  the  velocity  potential  resulting  from  'forced  motion'  with 
no  waves  present.  In  the  present  work  we  shall  consider  an 
initial-value  problem  for  ship  motions  with  forward  velocity 
and  show  how  Cummins'  (1962)  results  can  also  be  derived  from 
this  treatment  of  the  problem. 

We  shall  consider  the  motion  of  a  freely  floating  or  sub¬ 
merged  body  which  is  moving  with  a  constant  forward  speed  and 
oscillating  arbitrarily  (not  necessarily  periodically)  in  any 
of  the  six  degrees  of  freedom.  The  body  is  supposed  to  have 
zero  average  translational  and  angular  velocity  in  the  oscilla¬ 
tory  motion  in  some  appropriate  sense  of  the  word  'average'. 
Essentially  all  that  this  means  is  that  there  is  a  definite 
surface  moving  forward  at  a  given  speed  such  that  the  oscillatory 
motion  can  be  referred  to  this  surface  with  only  small  error. 

The  precise  meaning  will  be  explained  in  Chapter  II. 

The  position  and  velocity  of  both  body  and  the  free  surface 
are  assumed  to  be  knowr  at  some  fixed  instant  of  time  which  we 
shall  take  to  be  t  ■  0.  Besides,  we  allow  the  possible  presence 
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o £  waves  which  may  diffract  upon  tha  body  and  cause  it  to  move. 

The  Incident  waves  era  assumed  to  be  'known'.  Otherwise,  the 
nature  of  the  seaway  ie  eupposed  to  be  arbitrary.  Our  funda- 
isantei  assumption  Is  that  the  problem  can  be  linearized,  in 
the  aenae  that  tha  osclllationa  are  small  and  that  the  dis¬ 
turbance  of  the  free  surface  due  to  the  forward  motion  of  body 
la  also  email.  In  order  to  achieve  the  last  requirement  real¬ 
istically,  the  body  may  be,  for  example,  either  thin,  slender, 
flat,  or  deeply  submerged.  Aside  from  this,  w®  require  only 
that  the  form  of  the  body  be  'smooth';  no  symmetry  of  the  body 
is  assumed. 

As  usual,  we  shall  assume  the  fluid  to  be  incompressible 
and  inviscid  and  the  motion  of  the  fluid  to  be  irrotationai. 

The  analytical  method  which  is  used  in  the  present  work  was  first 
introduced  by  Volterra  (1934)  for  solution  of  certain  initial- 
value  problems  for  water  waves  and  was  later  extended  by 
Finkelstein  (1957).  Wehausen  (1965)  later  showed  how  such  a 
technique  can  be  modified  to  solve  a  class  of  problems  in  ship 
motions  and,  in  particular,  how  the  decomposition  of  the  velocity 
potential  of  Cunmins'  type  may  be  conveniently  made  by  this 
treatment.  The  present  work  is  an  extension  of  that  of  Wehausen 
(1965)  to  include  forward  motion. 

For  the  purpose  of  the  linearization  of  the  boundary  con¬ 
ditions  at  both  the  free  surface  and  the  hull- fluid  interface, 
we  keep  two  perturbation  parameters  in  mind:  £s  ,  measuring 
the  smallness  of  the  free-aurface  disturbance  caused  by  the 
forward  motion,  and  SH  ,  measuring  the  smallness  of  the 
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osclllatory  motion.  However,  in  order  to  have  a  development 
which  is  simultaneously  applicable  to  thin  or  slender  ships, 
and  to  deeply  submerged  bodies,  we  do  not  follow  the  traditiohal 
scheme  of  linearization  of  introducing  separate  boundary-value 
problems  for  each  of  the  order  £t  ,  etc.,  respec¬ 

tively.  The  boundary  condition  at  the  wetted  hull  is  linearized 
by  means  of  a  Taylor's  series  expansion  of  the  potential  function, 
in  which  all  the  terms  of  orders  £t  ,  £„  and  £*£*  are  kept. 

The  result  is  the  condition  which  is  also  called  the  Timman- 
Newman  boundary  condition.  A  somewhat  different  derivation  of 
this  boundary  condition  is  presented  in  Chapter  II  of  the  present 
work  so  that  additional  insight  into  the  nature  of  the  derivation 
may  be  gained.  As  to  the  free-surface  boundary  condition,  terms 
of  order  higher  than  £,£M  are  discarded,  so  that  the  traditional 
homogeneous  free-surface  boundary  condition  is  obtained.  For  the 
justification  and  limitation  of  the  applicability  of  such  a 
development,  we  refer  to  the  discussions  in  the  following  papers: 
Timman  and  Newman  (1962),  Newman  (1965),  and  Ogllvie  (1964). 

The. problem  is  first  formulated  for  the  general  case  of  un¬ 
steady  average  forward  speed  and  its  appropriate  Integral  equation 
for  the  velocity  potential  is  obtained.  Howevdr,  in  order  that 
the  decomposition  of  the  velocity  potential  for  the  'forced 
motion'  may  be  conveniently  made,  it  is  necessary  to  assume  the 
average  forward  speed  to  be  constant.  Cummins'  (1962)  develop¬ 
ment  for  the  case  of  constant  forward  speed  was  later  modified 
by  Ogilvie  (1965)  in  order  that  it  should  satisfy  the  Timman- 
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Newman  boundary  condition.  Hence  we  refer  often  to  Ogllvle's 
(1965)  work,  and  effort  la  also  made  to  preserve  the  same 
notation  whenever  it  is  convenient  to  do  so,  so  that  cross 
reference  between  the  two  may  be  easily  made. 

Uniqueness  is  not  established  either  for  the  boundary- 
value  problem  or  for  the  integral  equation  which  is  constructed. 
This  may  be  shown  in  both  senses  if  there  is  no  forward  motion. 

It  would  be  desirable  to  establish  this  in  the  present  case 
also.  Furthermore,  no  existence  theorems  have  been  established 
for  solutions  of  the  Integral  equations.  Thus  in  a  certain 
sense  the  work  is  purely  formal.  However,  if  one  is  willing 
to  concede  that  both  uniqueness  and  existence  should  be  provable, 
the  final  equations  show  the  proper  form  of  the  linearized 
equations  of  motion  and  the  nature  of  their  Ingredients.  In 
particular,  one  should  note  that  they  are  a  set  of  six  coupled 
lntegro-dlfferentlal  equations.  We  have  not  attempted  to  find 
any  solution  corresponding  to  a  special  geometry.  This  would 
be  a  reasonable  next  step. 


II.  Mathematical  Formulation 


Coordinate  systems .  It  will  be  convenient  to  consider 
simultaneously  three  right-handed  cartesian  coordinate  systems. 

Let  Oxfy  be  fixed  in  space  in  such  a  way  that  la 
in  the  direction  of  the  forward  motion  of  the  body;  0$  is 
directed  oppositely  to  the  force  of  gravity  and  the  (  J)- 

plane  coincides  with  the  undisturbed  free  surface.  The  coor- 
dinate  system  0  *  ^  J  will  be  taken  to  be  fixed  in  the  body  in 
such  a  way  that  when  the  body  is  at  rest,  the  axis  Of  i* 
directed  oppositely  to  the  force  of  gravity  with  the  center  of 
gravity  of  the  body  lying  on  the  line  of  ^ -axis ,  Ox  towards 
the  bow  and  0  j,  to  the  starboard;  and  when  the  body  is  at  rest 
the  (x,^)-plare  coincides  with  the  undisturbed  free  surface. 
Finally,  we  Introduce  the  coordinate  system  Oxjj,  ,  moving 
at  a  speed  equal  to  the  average  forward  speed  of  the  body  such 
that  when  the  body  is  at  rest  the  two  systems,  and  Ox'fJ. 

coincide  with  each  other.  Hence  the  two  systems,  and 

°*n  are  always  parallel  to  each  other;  0  x  coincides  with 
Ox  and  the  (  *  ,  j,)-plane  coincides  with  the  (A  ,  J  )-plane. 

In  particular,  if  we  assume  that  at  the  initial  instant  the 
two  systems  0x}£,  and  OX'Jf  coincide,  then  at  any  later 
instant  f  we  have 
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where  C(*>  la  Che  average  translatory  speed  of  Che  body. 

In  conformity  with  the  assumption  of  the  small  oscillatory 
notion  of  the  body,  we  shell  assume  that  the  displacements  of 
the  coordinate  system  from  are  small. 

Let  (*..  J.)  and  (.6,,  &tl  )  describe  the  linear 

and  angular  displacements  of  Ox'jJ  from  Oxyj.  ,  where 
*  1,  Thus  they  describe, 

respectively,  the  surging,  heaving,  sway,  rolling,  yawing 
and  pitching  motions  of  a  ship.  Note  that  at  any  instant  +  , 
the  position  of  the  origin  0  is  given  by  (  *i  •  1,  *  3r, )  in  the 
system. 

Suppose  that  &  ,  and  §a  a**  the  three  unit  coordinate 
vectors  of  the  frame  and  £  ,  £2  and  £are  those  of  the 

0?$J- frame.  If  P  is  a  point  in  the  body  with  coordinates 
(*,£,})  «nd  (X  >  ^  , £  )  when  referred  to  the  3Xft.  and  the 
-frames,  respectively,  then,  since  OP  *  QE  ~  0$  , 

we  have 

X  «I  +  f  1»  +£  2,  m  (x-x,)  e,  +  ■ 

Suppose  that  at  an  instant  t  the  body  frame  has  angular 

displacements  (  &,  ,  0Z,  $3  )  relative  to  the  Ox  f  y  -frame. 
Then,  without  assuming  smallness  of  the  angular  displacements, 
one  can  establish  the  following  transformation  between  the  two 
coordinate  systems: 
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X 

1 

X 

( K  0  o' 

COS$t  0  Sin®,' 

Cosa4  -Singj  0 

M, 

0 

o  cess, 

o  i  o 

S:"S*  COS® j  o 

> 

(2) 

lo  Sin®,  cotS, 

-Sin®,  0  tes StJ 

0  0  1 

A 

> 

. 

If  we  now  assume  that  &t  ,  and  are  all  small  enough 
so  that  one  may  assume  CosftZ  i,  s.-nfts  fy  ,  etc.,  then  (?) 
simplifies  to 


X  -  X,  =  X  4  0,  %  -  £3  y 

1- 1,  •  f  +*»*  -  fit  ► 

J-  *,  *  $+&•}-  &z* 


(3) 


or 


X-  X.,  =  X  + 


4.  Z,  J 


(3a) 


In  (3a)  we  have  used  the  usual  cartesian  tensor  notations  and 
the  repeated  indices  imply  summation.  If  we  now  solve  (3) 
for  X  ,  )  and  ^  and  assume  that  X,  ;  y,  and  have  the 
same  order  of  smallness  as  ,  $z  and  ^  ,  then  one  may 
discard  terms  of  the  higher  order  and  obtain 


■  ^ _ u 
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X  «  x-x,  -  0,|  4  fij  ft 

$  -  l-l  -$tx  +  t>>\,  ’ 
}  •  -  &>}  *  *1*, 


(4) 


or  X>i  *  ~  X;:  "  Zm  &j  t  (4a) 

Not*  that  the  unit  coordinate  vector  €>  >  for  example,  has 
the  coordinate  (X,  J  ,  J.  )  ■  (1,0,0)  in  the  -frame.  If 

we  let  J2/*X£>  *$£*$£,>  then  from  (4)  one  can  easily  find 
that  X  -  i  ,  $9-6}  *nd  $  »  <9,  ,  where  we  have  set  X,  =  Jf,*  },s.  0  , 

Hence  &  •  £  ■+(-03)£j  +  0,  £  ,  from  which  one  obtains  the 

following  result  which  will  be  useful  later: 

Sri*  1.  &•£  =  ,  *-1  •  &2  . 

Similarly  one  may  further  obtain  that  £i‘£>  *  1  , 

§i-£='0t>  Sri  3  4  and  §j'£i  ~  l  • 

Or  in  tensor  notation  we  write  this  result  as 
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Geometrlcal  description  of  the  ship.  The  surface  of  the 
body  will  be  given  at  all  Instants  In  the  body  coordina  te 
system  by  the  equation 

<6) 

By  the  transformation  (4)  this  function  can  be  written  as 


F(*-M)  -  F  •  F (*'!■}■*) . 

Thus  F  (*'$'%)*  0  and  F  (*' 0  .  respectively,  describe 

the  same  bpdy  surface  at  Its  instantaneous  position  in  the  body 
frame  Ox and  the  translating  frame  Ox ^  j,  .  Let  us  denote 
this  body  surface  by  S  .  Note  that  with  reference  to  the 
Ox  -frame  there  is  an  Imaginary  surface  S«  given  by 

S. F  (*'}'!)  -  0,  (7) 

which  is  stationary  with  respect  to  the  -frame  and  coin¬ 

cides  exactly  with  the  body  surface  S  when  the  latter  is  in 
its  undisturbed  position.  This  imaginary  surface  S<,  will  be 
called  the  "reference"  surface. 

Description  of  the  sea.  The  form  of  the  free-surface  Ox^J, 
will  be  described  either  in  the  space  reference  frame 
by  the  equation 

}  -  ?  (*■!■*)  , 


(9) 
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or  In  eh*  translating  reference  frame  0*}}  by  the  equation 

?  =  Y (*■}'*) ,  <»•> 


where  we  let 

y  (X'Z't)  -  ?(*+£*' *.t)  *  y  (*’},+) . 

Potential  functions  and  their  preliminary  decompositions. 

As  usual  we  shall  assume  the  fluid  to  be  heavy.  Incompressible, 
and  lnvlscld,  and  the  flow  to  be  Irrotatlonal,  so  that  a 
potential  function  may  be  defined.  Let  (%’  V1‘  1)  be  the 
potential  function  such  that  its  gradient  equals  the  velocity 
vector  of  a  fluid  particle  with  respect  to  the  space  reference 
frame  .  This  velocity  will  be  referred  to  as  'the 

absolute'  velocity  of  a  fluid  particle.  Me  shall  also  write 

The  relative  velocity  of  a  fluid  particle  with  respect  to  the 
translating  reference  frame  then  should  be  given  by 

V  *  =  (-c  +  fc,  (fv  (ft)  .  (9) 

Note  that  (fa  s  (fK  ,  4}  -  ^  etc .,  but  $  *  “  c  % 

V  (f  still  gives  the  absolute  velocity  but  is  expressed  in 
terms  of  the  variables  of  the  system. 
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In  the  subsequent  development  we  shell  suppose  thet  the 
resultant  fluid  motion  of  our  problem  is  composed  of  the 
superposition  of  two  parts:  1)  the  disturbance  due  to  the 
translation  of  the  body  fixed  in  its  undisturbed  position  with 
the  forward  speed  C(+)  into  otherwise  undisturbed  fluid  region, 
and  2)  the  fluid  motion  due  to  the  oscillatory  motion  of  the 
body  and  the  oncoming  waves  Hence  we  shall  write 

(*•?•},*)  =  %  +  ft  (*•!'}.+)  ,  (10) 

where  %  and  (f  ,  respectively,  represent  the  fluid  motions 
due  to  the  first  and  the  second  parts  mentioned  above. 

The  linearized  kinematic  and  dynamical  boundary  conditions 
on  the  free  surface.  A  systematic  linearization  of  the  math¬ 
ematical  expressions  for  the  fact  that  the  free  surface  is  a 
material  surface  and  the  assumption  that  the  pressure  every¬ 
where  on  the  free-surface  is  constant  will  lead,  respectively, 
to  the  following  conditions; 

^  ■=  o,  (U) 

and 

(I2) 

They  are,  respectively,  the  kinematic  and  dynamical  boundary 
conditions  on  the  free-surface  and  are  to  be  satisfied  at  the 
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undisturbed  free  surface  ^  ■  0  .  Elimination 

of  Y  from 

(11) 

and 

(12)  gives 

on  £  =  0. 

(13) 

The 

relationships  J  >J,t )  *  , 

Y(*Zfi)m 

YMh 

<h 

*  etc.  and  =  tyt~c  fy,  give  us  easily  the  counterparts 

of 

(11),  (12),  and  (13)  in  the  0*11  system  as 

i  follows: 

( ?rci*)Y( w)  - 

=  0, 

(14) 

f  *  (if  ~c$<)  9 

-  0, 

(15) 

and 

(h-chi  +  }<j>}=c 

>  > 

(16) 

whlch  are  to  be  satisfied  on  the  undisturbed  free  surface  J*0  . 

Linearization  of  the  boundary  condition  at  the  hull-fluid 
interface .  The  boundary  condition  at  the  hull-fluid  interface, 
like  the  case  of  the  free-surface  boundary  condition,  can  be 
formulated  either  in  the  space  reference  frame  or  in 

the  translating  reference  frame  Oxty.  It  will  be  convenient 
in  the  future  if  we  work  here  with  the  translating  reference 
frame  Oxf}  . 

Let  us  first  work  out  an  expression  for  the  velocity  of 
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the  ship  hull  relative  Co  Che  crane  lacing  reference  frame  0*jf 
according  Co  Che  transformation  (3)  and  (5)  valid  for  Che 
small-oscillations  approximation.  Let  (J  be  a  typical  point 
on  Che  surface  of  Che  body  with  coordinates  (  })  in  Che 

body  reference  frame  Ox  and  in  the  frame. 

Then  the  vectors  X?  »  X  £<  +  $  h  $  !=j  =  0  Q  end  X  = 

X  &  4  f  m  OjH  are  the  position  vectors  of  Q  in 

end  Oxyj,  respectively.  It  is  easily  verified  from  (5) 

that 

h  -  e,  -t  e,  .  a?) 

Hence  we  may  write 


x'  .  S  £  +  7  £,  4 1  £ 

&■+$  e,  +  + fit )e,  4 

4  h  . 

Let  us  define  2  "  X  ± /  4  4  £  £}  ;  then  we  have 

X'  =  X  4  *  x  >  (18) 

A 

Obviously  X  defined  above  is  not  a  position  vector  of  the 
point  Q  with  reference  to  any  of  the  two  reference  frames 
considered  here.  We  may  now  write 

X  =  ,  X,  4.£  4  &3)X  £  .  (19) 
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The  formula  (19)  gives  immediately  the  velocity  of  the 
point  4  relative  to  the  0*fJ,-frame  as  follows: 


V (<?;*)  *  x(4;t)  *  +  ,  <20> 

a 

where  we  take  £  s  0  since  the  body  is  to  be  assumed  rigid. 

If -the  vector  ^  in  (20)  is  resolved  along  the  Ox]}- frame 
we  may  then  write 


V,  (*  *1 


X«  <+>  +  Sj  Xi 


(20a) 


Let  J2  be  the  unic  normal  vector  to  the  hull  surface  point¬ 
ing  into  the  body.  Let  JJ  be  resolved  along  both  reference 
frames  so  that 


n 


=  i  &, . 

r,  t 


(21) 


Suppose  that  the  function  p  given  in  (6)  which  describes  the 
body  surface  is  so  chosen  that  the  inward  normal  is  given  by 
the  formula 


A 


(22) 


where  ft.  represents  the  partial  derivative  of  the  function 
F  with  respect  to  its  i-th  variable.  Then  from  (21)  and  (5) 


we  obtain  easily  that 
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*  \?>i  +  £>*-•  h  Pi  J  /j?>*K  .  (23) 


The  boundary  condition  on  the  body  is 

n  •  v(-cx  +  ({  )  =  a-]/ 

i.e.,  the  normal  component  of  the  fluid  velocity  at  a  point 
<3  on  the  body  surface  equals  that  of  the  surface.  In  component 
form  this  can  be  written  as 

-  .  (24) 


Note  that  in  (24)  variables  of  both  coordinate  systems 
are  Involved.  Transformation  (3)  can  now  be  used  in  writing 

^(*•1'}'+)  *  4(*  +  kl+hS-9t$'  St*  ’t) 

Since  the  quantities  X,;  and  Qx  are  small  compared  to  the 

A 

quantities  Xt  ,  it  is  reasonable  to  assume  that  the  following 
Taylor' 8  series  expansion  for  the  potential  function  fy.  is 


possible: 


16- 


&(**'*’*)•  fyA* •Ht)+[Xu+Zim*0„X"lbu+-’  ‘  .  (25) 

This  is  now  to  be  substituted  back  into  (24).  One  may  at  this 
point  Introduce  several  perturbation  parameters,  for  example, 
•ay,  £t  measuring  the  smallness  of  the  free-surface  disturbance 
due  to  forward  motion  and  measuring  the  smallness  of  the 
oscillatory  motion.  Then  by  following  usual  scheme  of  lineari¬ 
sation  one  can  deduce  from  (25)  and  (24)  linearized  conditions 
for  separate  boundary- value  problems  of  the  orders  £,  ,  and 
£,£„  ,  etc.,  respectively.  However,  the  goal  of  the  present 
work  will  be  better  served  by  following  a  slightly  different 
approach.  We  shall  not  be  so  specific  about  the  Introduction 
of  the  perturbation  parameters  but  rather  shall  discard  what¬ 
ever  terms  of  order  Qj  (9*  or  Kj  appear  in  (24) ,  for  they 
are  clearly  terms  of  the  order  .  (Terms  of  order  £r £ m 
will  Always  be  retained.)  If  this  is  done  then  we  have,  from 
(24)  and  (25), 

K  (»■}■}>  &<*■  %U)  *  %  i  ■*  £, 

/\ 

where  the  expression  (22)  has  been  used  for  /?„•  .  If  we  now 
drop  all  the  circumflexes  over  the  variables  in  the  last  ex¬ 
pression,  we  have 
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F*(*W)  tyu  (*•)•}'*)  Qj  F'kfo  ^  FA  K(  +ltm6mK]  fysi  - 
=  PiAx.tVlkt  i^iAfyxA].  (26) 


Obviously,  (26)  is  now  a  condition  for  the  potential  function 
(f  to  be  satisfied  at  a  point  on  the  surface  defined  by  the 

A 

equation  f  0  which  if  precisely  the  imaginary  sur- 

face  defined  by  (7). 

Although  certain  physical  interpretations  of  the  implication 
of  (26)  are  possible,  we  shall  not  do  this  but  rather  refer  to 
the  original  paper  of  Tinman  and  Newman  (1962)  and  to  the  dis¬ 
cussion  on  p.  39  of  Ogilvie  (1964).  However,  we  should  like 
to  remark  here  that  the  fact  that  (26)  is  to  be  satisfied  at 
an  Imaginary  boundary  S,  comes  out  naturally  as  the  result 
of  linearization  and  the  way  (f>  is  expanded  into  Taylor's 
series  in  (25).  The  question  of  whether  S.  represents  the 
mean  position  of  the  oscillating  surface  S  or  not  is  immaterial. 
In  fact,  S.  here  will  seldom  be  the  mean  position  of  5  , 
since  we  are  considering  an  arbitrary  oscillatory  motion  of 
a  body  in  forward  motion. 

Let  us  now  put  (f  *  (ft  ■+  (f>t  in  (26)  and  assume  that 
(f>  ,  and  its  derivatives  have  the  same  order  of  magnitude 

as  Oj  or  X,;  .  Hence,  terms  like  Ejjt  6>j  F.i  and 

mnfimXn)  ^  are  to  be  discarded.  Then  from  (2b) 
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we  have 

Fjnrtt/x-U'f)+t  (*'>■**>! -+  + 

+  h  [*./  4  hm  K  'C&1&  +  l'ii  *i  ?+  1 

=  L<*‘U)[x,i 4  &#*4**].  (27) 

Suppose  that  there  Is  no  oscillatory  motion  (i.e.,  we  set  ?h~°)> 
so  that  XtS  s  0,  #.«0,  and  ft  •  o  ,  ft  =  ft*  .  Then  from  the  last 
condition  we  have 

F-sM'Dfaih'W.*)  -c4  &  -*• 

Hence  ftaj  ft  .  -  C  fl0l  on  S,  :  F  (xrJ ■}■ )  -  0  ,  where  we  have 

written  HtjQt)  s  F.i  ($•)//&£*>  for  £  S#  .  Or  we  write 

(28) 

With  this  result  another  condition  for  ft  can  now  be  obtained 
from  (27).  Thus  we  have,  from  (27)  and  (28), 


L  (*-n)  i;  (»bb*) 4 Si  F,t  ft"; 4  L  [*,, 4  u . &...K]  fttil 

~  c  £*j*  4 
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Let  us  again  write  flaji  (*■}>})  for  F-~  (*' and 

rearrange  terms  in  the  last  expression.  Then  we  have 


-^0*  <29> 


For  convenience  let  us  introduce  a  vector 


defined  by 


A  =  X  -f  .  (30) 

Or  in  component  form,  ^k  •  From  (30)  we 

have ,  then , 

If  *  x,  +  (6,4,4)**  >  (3D 

which  is  a  vector  evaluated  at  a  point  Q,  *  (*'|f<$.)on  S, . 

Note  that  the  difference  between  the  vector  given  by  (31)  and 
that  of  (20)  is  of  the  second  order.  Hence  (31)  may  be  re¬ 
garded  as  the  first-order  approximation  of  the  velocity  vector 
of  a  point  Q  on  the  actual  body  surface. 

Condition  (29)  can  be  put  into  a  compact  form  if  the 
following  vector  identity  is  used: 


v*L6*&]  -  U-V)&-  (a-v)±*AV-A-±r-A  ■  w 
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Let  X  »  V  (-CX+  fa)  ■+  y  (f>0  and  A  be  as 

defined  by  (30);  then 

1)  7*2  *  V-V(-c*+b)  =  V'C-cx+faoo, 

2)  V’A  =  P*  /*• 4  ^ £ ) 

*  ty  4;  *  ty  =  o, 

3)  te-pj-d  *  f[t 7(’C*+ft)]-r)l*/  +  (&>b.4)*.X } 

*  )  [  *v  4  £->*  **  7,4 

=  (’cfmht)L*vJL0jl 
yA  4*  4  ^  )  ; 

-  ^*7j  A  =  -  [[A,  +  6M,(  63)xX]-v}v(-^ifa) 
-HytqXtjku  . 


4) 
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Note  that  3)  and  4)  are  precisely  the  same  as  the  second  and  the 
third  term  in  the  right-hand  side  of  the  equality  in  (29).  Hence 
with  the  use  of  the  identity  (32),  (29)  can  now  be  written  as 

£.•  7#  ®  +  /7.*7>M  *V(-c*+to)}'  03) 

As  was  remarked  before,  to  the  first  order  approximation  we 
have  9^/6'/  =  {/  ,  the  velocity  of  a  point  Q  on  the 

instantaneous  position  of  the  hull  surface  relative  to  the 
OAjj,- frame.  We  shall  then  rewrite  (33)  as 

(f  U:t)l  =  y„  (i;f) 

ii  n  Is, 

(34) 

where  U  -  V  (~C  *  + <jf>.  (*,},}, 4))  .  and 

x  —  ( Xi  y, $ )  1  So . 

Note  that  (34)  is  precisely  the  condition  originally  de¬ 
rived  by  Timman  and  Newman  (1962).  It  might  be  interesting  to 
note  that  in  (34)  an  additional  second  term  is  needed  in  a 
condition  derived  from  the  kinematic  boundary  condition  on  the 
hull  surface  S  .  However,  it  should  not  be  surprising  that 
this  should  be  the  case,  since  S„  ,  being  an  imaginary  boundary, 
is  not  a  material  surface.  Therefore,  even  though  *  V* 

is  the  appropriate  condition  for  a  material  surface,  it  need 
not  hold  on  S»  .  and,  in  fact,  when  a  body  is  translating  as 
well  as  oscillating,  it  does  not  hold. 

Finally,  if  the  condition  (34)  is  specialized  to  thin-, 


mm . m  ■■■■■  i  '  ■■■■ . ea 
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flat-  or  alandcr-ahlp  approximations,  it  is  found  that  (34) 
in  fact  rapraaants  a  combination  of  several  boundary  conditions 
for  separata  boundary- value  problems.  Moreover,  in  this  case, 
the  additional  sacond  term  contains  terms  of  higher  order  in 
perturbation  parameters.  Such  terms  should  be  discarded  if  a 
strictly  firat-order  expression  of  the  boundary  condition  is 
required.  Nevertheless,  the  present  form  of  the  boundary 
condition  in  (34)  willbe  used  in  the  present  work  because 
there  seems  to  be  no  harm  in  keeping  some  of  these  higher- 
order  terms,  provided  it  is  understood  that  the  accuracy  is 
good  only  to  the  firat-order  approximat ion. 
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III.  Solution  of  the  Initial-Value  Problem  by  the  Method  of 

Green's  Function. 

The  initial-value  problem.  In  the  subsequent  development 
we  shall  need  only  the  translating  reference  frame  Ox'fy. 
Henceforth,  the  surface  5.  •*  0  will  be  restricted 

to  the  part  of  the  surface  which  is  below  ^  =  o,  i.e.,  $o 
coincides  with  the  wetted  part  of  the  hull  surface  when  the 
ship  is  in  its  undisturbed  position.  We  shall  be  looking  for 
an  unsteady  velocity  potential  ^  (*'),}•  i)  satisfying  the 
following  equations  and  boundary  conditions: 


&  +  f *°  ■  }<0‘  <35) 

<3‘> 

on  the  undisturbed  free  surface  ^*0; 

(37> 

S»  *  7 

where  X  *  (x-  in  S,  ,  and  V»  (&  ’i)  -  V(~£x  4 
already  prescribed; 


G.  |8  •  0 


(38) 
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where  B  is  Che  bottom.  If  the  fluid  is  infinitely  deep,  the 
lest  condition  is  repleced  by 

(p  )  =  0  .  (38a) 

In  addition,  ,  (f^  end  their  first  derivatives  are  assumed 
to  be  uniformly  bounded  at  °0  . 

The  initial  position  and  velocity  of  the  ship  are  assumed 
known.  Moreover,  at  the  initial  Instant  both  Yi  and  Yit  on  the 
free  surface  are  prescribed: 


Y (H>o)  --  f,  (*■})  t 

Y*  (**}•<»•  A  (x.}) 


(39) 


where  -f>  and  ~fi  are  given  functions  of  X  and  y  . 

From  (39)  and  (14)  we  see  that  is  also 

determined.  Thus  at  i*0,  fa  is  given  on  all  boundaries  and 
is  bounded  at  infinity  and  therefore,  °)  can  t>e 

obtained  as  the  solution  to  a  Neumann  problem.  Ue  also  know 
from  (15)  that  Y>  (*'}”*)  ,  the  free-surface  elevation  due  to 
the  fluid  motions  associated  with  ,  is  given  in  linearized 
theory  by 


Yi  (*'}•*)  -  ~J  (jT'C  If)  (ft  (x  .o,yt)  . 


(40) 
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Hence,  from  the  initial  values  given  in  (39),  tyt 
and  $ tt  can  also  be  found. 

The  time-dependent  Green's  function  and  use  of  Green's 
theorem.  The  Green's  function  G  (*'  </<i  <  f- '/>?;  4  )  in  question 
is  required  to  be  a  solution  of  Laplace's  equation  in  the 
variables  7.  £  with  a  singularity  like  !/\r  * 

■to 

*  [ *t  (*'}•}■)  but  otherwise  harmonic 
in  ? ,  'j,  f  in  the  region  ,  >)  *  0  .  Thus 

G  (x,nn>'/>  ?;/;  =  */r  i  H  1.J,  m)  (41) 

where  f- '/,  ?;  f )  is  harmonic  everywhere  in  the  domain 

of  definition.  We  shall  suppose  that  G  is  a  symmetric  function 
of  i  satisfying  the  following  boundary  and  initial  conditions : 


(hmenhf  ■*  }G?  =  0  . 

G-v  (-  ‘  1 1  ~t)  |b  *  0  or  Gj  -  0  , 

G(z>  0)  =  0  ,  GtU;  ?'Pi;o)  =  0  , 
G  *  & (-R  )  >  G/i  -  O'  (R3)  as 


(42a) 

(42b) 

(42c) 


(42  d) 


26 


where  R  *  i  and 

sT  *  Hi  (l)  • 

A  mat  hod  of  construction  of  such  Green's  function  can  be  found 
aithar  in  Stoker  [1957,  pp.  188-191]  or  in  Wehausen  and 
Laltone  [1960,  pp.  491-495].  It  is  a  property  of  this  Green's 
function  that 

G(miit)  -  G(ti4Li*)  (43) 

where  (*’?.£)  and  1  ■  ft'/*/)  ia  a  point  on  the  boundary 
of  the  fluid  region  considered. 

With  Green's  function  described  above  we  may  proceed  now 
to  set  up  an  integral  equation  for  the  function  ^  .  To  achieve 
this  goal  we  start,  in  the  usual  fashion,  with  applying  Green's 
theorem  to  the  Green's  function  Gf  and  to  in  the  fluid 
region  bounded  by  the  undisturbed  free  surface  F  ,  the  re¬ 
ference  surface  $«  ,  the  bottom  3  (if  any) ,  and  a  large 
sphere  H-  of  radius  <1  centered  at  the  origin  of  the  trans¬ 
lating  reference  frame.  Note  that  only  parts  of  f  ,  B  and 
II  will  serve  as  bounding  surfaces  and  we  shall  call  these 
parts  F  t  B '  *««*  -ft/.  Then 

fyfteii)  *  [j  /  ;  I; *-rJ  4^(1  ;*)-  $v}  </s,(44) 
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where  the  normal  vector  Is  taken  to  be  exterior  to  the  fluid 
region  considered.  The  surface  integral  over  J)1  vanishes  as 
the  sphere  JL  extends  to  infinity  because  of  the  boundedness 
of  and  fyf  and  the  behavior  of  G  at  ^0  .  The  integral 
over  ft'  also  vanishes  since  both  and  G)v  are  zero  on  [ft  . 
After  letting  Q  -*  ,  (44)  then  becomea 

4/f  1 *)  •  4ir  H I Cl (JL '  (t  °’**)- titty] 

*  4*  II I  ;1 ;  ^  fev  (li  V~tiltG~)<Js.  (45) 


Interchanging  i  with  T  in  (45)  and  observing  the  fact 
that  G  (&}t  !  4-t)~  G  (£:1;  T~t)  ,  let  us  now  Integrate 
both  sides  of  (45)  with  respect  to  T  from  0  to  f  .  Then 
we  have 

4lt  ({,(*;+)'  4*  ti  U;o) 

=  /  Jr  jj j £(&>%. o,t;4-T)  </>  (w.t)~  J 

*  F  ' 

♦A// £(*•  J  '■  }  <4 

*  *  J>,  • 


(46) 
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We  may  proceed  further  with  (46)  by  making  use  of  the 
boundary  condltiona  aa  follows.  From  (42a)  we  have  now,  on 
the  undisturbed  free-surface  F  , 

t)  *  -  -f[  (&  ~c(*‘T)7f)G(i;to,x;4-Tj\  (47) 

where 

h  £(£‘1 ;  *-t)  2  <7/ - : ? '  * •  Pr  G  ^ 


is  the  derivative  of  £7  with  respect  to  Its  seventh  variable. 
With  further  use  of  the  free-surface  boundary  condition  for 
the  potential  function  (ft  ,  it  is  not  difficult  to  verify  that 
the  following  Identity  holds  for  1  *  ('  }-  i  )  on  f  . 

*  7f  f^L^y + y> (- /~r)j f  ■ 

~i  Tf  I  *  t-rfcfo+W-riJ&y 6  (* >*-*)) 

+  lk<W'-r)lCM(*6)]} 

4  ■j’  i  J^r  |^t)  -f^Cfr)  ■>  Cf/-r)J  jiipQ)  ■* 

4 [*(Tl-c <*-#](}(,)+ cMLca-v-ch,]^  fcjJ  f 


V  (43) 
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+  j  ^  n)\z  CM ($-lG)  +  c(r>[c(T)  +  2  c(4-t)]£(j>6)- 

-c(t)(DG))1 

C(4-T ) 


where 


J)Y,( rt.f)  *■ 

<Rd  DGU;  1;4-t)s  -  c  (+-•*)  G(*,l ,  *-tj  . 

With  the  use  of  the  Identity  (48),  the  Integral  over  the 
undisturbed  free  surface  F  1ft.  (46)  now.  becomes 

if  -  -//Is [Mt to)]  +  ]}Jt4{ 

F 
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where 


F 

-  [3  c (T)  ic a  -  T)]  ft  (M)  i  [c (T>- C  tt-T )](D$y 
+  C(T)[c«-t)-c( T)]7f(PG)j 
4  ^  l  ^T/  2.C(r)  jfQ  V-t)  + 

-t  cm [c  (-cm c a- r>j  j^(DG)-  c(t) (2>6]j  ^ 

jt  </, t  (l  o,  J, T)  j  &T)  D G,  (* )  1. 0,1 ;  i-T)  }  </  f 

+  J ^  <V l  T){Cl ^‘T)Tf  $  (* 1 1  *' 'l ! .  > 

F  1 

DGUilM'.o)  H  (£-<(0>7f)  G(£;t°*t;o)' 


In  deriving  (49)  we  have  made  use  of  the  condition  (42c), 
namely,  £  (X  ,  ?,  O,  $  ;  o)  =  0  ,  (&  !  I  f.f;  0 )*  0  .  it  follows 

from  this  condition  that 

}G(*;iM;0)  *  Gt(z,  ?, J, *; o)  (*  ; t °<{, »)  =  ^  . 

Note  also  that  in  obtaining  (49),  Gauss's  theorem  for  the  plane 
region  has  been  used.  Since  the  undisturbed  free  surface  F 
may  be  considered  as  being  bounded  inwardly  by  P„  ,  the  contour 
around  the  intersection  of  £T0  with  f  ,  and  outwardly  by  /\> 
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at  infinity,  one  of  the  integrals,  for  example,  becomes 

-  i  ^(ro.i.r) [e-ftj  2)6 ]Jr 

I  4 

where  f*  i»  oriented  couterclockwise,  and  the  integral  around 
r »  vanishes  because  of  the  behavior  of  £  and  the  boundedness 
of  ^  at  infinity.  The  other  contour  integrals  appearing  in 
(49)  are  obtained  in  a  similar  fashion. 

We  shall  next  proceed  with  the  Integral  over  S.  in  (46) 
by  making  use  of  the  boundary  condition  of  the  function  ^ 
on  St  •  Let  us  first  Introduce  the  following  notations.  Let 


£/<  ,  -#*/.  2,  3, 

,  *m  4’s> 


(50a) 


Hi 


n.  (s)  •  € k  ,  i  *  1,2,  3, 

< 


i  =  4.S’,L, 


(50b) 
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!!»(*)•  j  V*  [tk*  ¥.(*;*)]}>  k-t't‘3. 


(50c) 


n.  (*)•  I  V*[  (c<.3  *£)  *  y.(*'+>]},  i  -4>sr,  U, 


where 


Vo  (x;1)  38  V  (-CK  -i(f0  (x it)) 


<  °<x  *  X  ,  <*,  *  },  ;  °U  -  *  °<s-~Sz,  °<i~  &3  .  <50d) 

Then  the  vector  A  defined  in  (30)  can  be  written  as 

A(l,i)  =  X,(t)  -t  (ffi,.*),.*,  )  X  X 

'  ^  M  (*)  •  (51) 

Thus  (37),  the  boundary  condition  for  on  S',  ,  can  now  be 
written  in  the  following  form: 

fyh(*u)l  -  X  i  t  Lfat)  tty*) .  (52) 

I  s#  *»»  * 


Taking  the  time  derivative  of  (52),  we  have 


-33- 


=  i  itdt) 

"Hf  s , 

( 


where 

ij(-i  4  )  =  ^  (-  '  * 

The  condition  (53)  is  now  to  be  substituted  into  the1  integrand 
of  the  integral  over  S*  in  (46).  If  this  is  done,  we  have  then 


*  jt  jj it  /'  t-r)  [ %(!)  M  -f  Z  (t(t;r)  <*^(t)  + 

h(i;-')*<i(T)  Ms 

-  If  ft  bn)  ds 

s. . 


>  (54) 

i 


JJ  ^>°)  £*(&■!£;  -t)  ds  -  Jj <f>  (t)  f)  Gitteili  °)  ds 

S.  s. 


J 


Note  that  (54)  end  (49)  together  equal  the  right-hand  side  of 
(46). 

Define  the  operator  as  follows: 

V 

J  {fy  }  (X;+)  2  An  h(x-,t)  ■+  jj  ty,  (ti +)£„(»:  Vo)  Js  * 

J  s. 


+  J  Jr  Jjty,  (l;r)  (Hill  li-T)  Js 


S. 


I 

7 


j fj tyi(1,0't‘r)^\C(T> 4~t)~ 

F 

-[3CMi  c(i-r)]^(DQ)  ■+  [  C(r)-  c(i-T)]  (!>£)+ 
+  C(t)[C(4-t) -c  (r) ]  ||  (T>6  )J 

/ tjrJj^(W'r)  ItAt;  /-t)J  J-fJs 


i)\l  C(r)  D%  I  i-r)* 

+  C(r)  [c  (T)-tlC  (i-T)]  i  ($G)  -  c  (t)(JH 5)}  Ji 

ll*Jr  <f  {g%(to.lT)j^(r)J)^U;l,o.ti*^  Jt  t 


(55) 


With  this  definition  of  d[  ,  the  following  integral 
equation  for  (ff  can  be  obtained  from  (46),  (49)  and  (54): 
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— 

*{<r,  }(»*>  =  47T  fy(x;o)  jjfy  (t;o)£J&;±;i)  ds 

s. 

Y.frGfri'O-r.-vljJfJf 

F 

4  f  [c(0HCCt)]  II  (ft(i 0, j, 0)  6  (£  i  * ^  i ;  o}  M  ►  (56) 

-j-C(o)j)  j} 

'  r§ 

i  J  Jr  nG (&’>  i  ;V- t)  I X  n,t(i)%t(-r)  + 

s.  ^  . 

+1  4ft;-'-,)  V7'  *  jr,  lJs. 


So  far  our  development  has  been  .perfectly  general  In  the 
sense  that  no  restriction  Is  Imposed  upon  the  average  forward 
speed.  Henceforth  we  shall,  however,  assume  that  c  “  const., 
l.e.,  the  average  forward  speed  of  the  ship  Is  a  steady  one. 

As  a  consequence  of  this  assumption  we  have  now  , 

^  =  »  *nd  ¥•  -  ,  l.e.,  they  become  Independent 

of  time.  Thus  the  only  time-dependent  functions  are  the  Green's 
function  £}(&i  £  ;i)  ,  the  unsteady  velocity  potential 
and  the  various  displacements  #4  (i)  .  The  integral  operator 
defined  in  (55)  and  the  integral  equation  (56)  now  become, 


respectively, 


o fib  }(*■>+)  =  (f,  (*;t)  +  U 

iJo  Js 

s. 

-  j  [ rJJ(f,,(r°'  r •  V  ^  l D2G< (*  t  o.t, ;  1 -t)  - 

F  -4&&G)  )Jt<h  j 

-cj'  drpj  )J$f  G(&  < 

r*' 

e  r,  j  ) 
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and 


JlfyjOni)  =  4irfy(jt;o)+Jjk(VO)&*(*:i;t)Js 

s. 

‘  //{&  [W(t  v;o)]+  Y[2G(»; 

^jlji  ft  **;«)/£  W*f  *•  *;  f ;  ]  it  ir 

"  ^h-O'l.O)  j}G(x-.1.o, T;  o)  ^ 

ij'  <tr  jjGUil-J-r)^  Kofi) -t 

*f  ^  faJ  ^(rjj  Js 


(58) 


We  see  that  integrals  over  the  undisturbed  free-surface 
/■  are  involved  in  the  expressions  on  both  sides  of  the 
equality  in  (58).  The  behavior  of  Green's  function  G  at 
infinity  of  our  problem  makes  these  integrals  converge  even 
though  the  domain  of  integration  p  actually  extends  to  in¬ 
finity.  In  case  of  zero  average  forward  speed,  i.e.,  c  ■  0, 
all  the  integrals  over  F  and  around  Po  in  the  definition  of 
the  integral  operator  disappear.  Hence  equation  (58) 
in  this  particular  case  can  be  made  to  yield  an  integral 
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equation  for  a  function  defined  only  on  the  reference  surface 
S,  by  letting  the  singular  point  *=  (*'}'}')  converge  to  a 
point  of  the  surface  S.  .  However,  with  non-zero  average 
forward  speed,  we  shall  inevitably  deal  with  an  integral 
operator  involving  both  surfaces  S#  and  f  .  Therefore,  a 
considerably  greater  difficulty  should  be  anticipated  in 
solving  an  integral  equation  obtained  from  (58). 

The  following  identity  may  be  established  from  (57): 


+  lj  ft  (Vo)  Js 

s. 

'  J H k  (*'  r >o)  %  l  ;  f' *' r ; 

F  -  4^(J>G)  J  JrJs 


r  (59) 


r i » 


From  (59)  and  (56)  we  find  that  must  satisfy 


-39- 


=  ~  IJ  i  £*[%(*' t*°)l 

F  -+V  [^.DGUiW:i)]]dtJK 
+  J  1/ fa  (*•*•*'  f 1 6  (& ;  t  o.  ?,  i)  ~  z&  (VG )}  M  4 
+  c  H  ^  (t  *-T.  o)  G(£:to,  r;+)}  Jf  Vr 

--j  <ft  ty,  (I'O.x.o)  \T>G)(&:'t'0‘V  >'?)J 

fo 

<jcf>  4lfr*T.o){nfrt«r;*)+SC% 0>6)- 

-&&£) }  ^X 

+  I/g  (X;i  >o)\i  n>/t) ifi)  «i(v }  Js 
*1  II ' 1 '  *r,  n9fi.)^(r)  +2 

S.  *“ 


It  would  be  desirable  to  establish  the  uniqueness  of 
solution  of  the  integral  equations  (58)  and  (60)  and  of 
equations  to  appear  later  which  are  of  the  same  form,  namely, 

JlU'f- 

That  is,  we  hsould  like  to  prove  that  implies 

that  Cf  -  o  .  Although  this  can  be  established  for  c  •  0, 
we  have  not  been  able  to  prove  it  for  c  i1  0.  However,  it 
seems  very  likely  that  this  is  so  and  we  shall  assume  hence- 
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forth  that  it  can  be  established. 

From  the  physical  situation  we  are  considering  here  and 
from  the  linearity  of  the  problem,  it  seems  to  be  clear  that 
the  unsteady  part  of  the  velocity  potential  has  the  following 
constitution: 

^  ■>  i, 

where  ^  ,  and  represent,  respectively,  fluid 

motions  due  to  a)  forced  oscillation  of  the  body,  b)  the 
incoming  waves  and  c)  the  diffracted  waves.  In  order  to 
consider  a  general  situation,  we  do  not  wish  to  assume  that 
the  body  starts  oscillating  from  a  state  of  rest  relative  to 
0x?pat  the  initial  Instant  4*0  >  on  the  other  hand,  it  is 
desirable  to  do  so  for  the  convenience  of  the  type  of  de¬ 
composition  which  we  shall  consider  in  the  next  chapter.  As 
a  possible  approach  to  solving  such  a  dilemma,  we  shall  assume 
that  the  velocity  potential  (fF  may  be  further  divided  as 
follows: 


-  %o  t  , 

where  tyFo  describes  the  fluid  motion  which  would  take  place 
as  a  result  of  only  the  given  motion  of  the  body  at  the  initial 
instant  and  fo/  represents  the  fluid  motion  due  to  the 
oscillation  of  a  body  which  has  started  from  a  state  of  rest 
and  has  achieved  the  given  Initial  motion  of  the  body  Instantly 
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at  /  —  O  .  Hence  tyFo  and  ,  satisfy,  respectively,  the 
following  boundary  conditions : 

%on(*;o)lsr  t(*;o)lSo.  KoJ*:t)ls0  =  °  for  iyo] 

and 

For  convenience,  let  us  henceforth  write 

#  -  fe,  +  6, , 

where  =  tyo  • 

The  two  functions  Cfn  and  (f^  for  i  >  0  are  then  defined, 
respectively,  by  the  following  equations: 

i  I  &/ )  :i>  -  lJTfl </>Jvr)  Js 

So 

=  ]  Jr JjG(x;1;S-T)\ij>.k(-i)«t(r)+ 

■f  Z  iA(t)  <kA( r )  J  Js  t 


(62) 
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and 


4-w^(jao)  *  jl(f,(r.o)G^(S:l:i)  h 

s4 

-JJjG[2)Kato)]-tY[3)G(x;  1,  o.xs  t)] }  *  h 
F  l 


(63) 


j>  ty/tO'X.o) 


fyi  *  and  ^  ,  respectively,  satisfy  the  following  initial 
and  boundary  conditions : 


(*  -,o)=0  ,  (fnt(x.o,^)  0)  =  O  , 

tyfin  (*  ,0)lso  ~  0  ' 

?,<»  M/s.  -  kJz-'ttl. 

=  i  LU!0L(h 

*=>  &*>  * 


(64a) 

(64b) 


(64c) 


for  i  >0 
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9 «,(*■••  o)  =  (f,  (*  'O)  . 


(65a) 


9*t(*-°'hO)  -  9lt(*'0'^°)  - 


(65b) 


*1  ItotttXktCo)  +1  l(&)\(°),  (65c) 
#•/  *•' 


C-  '  ^ 


=  0  for  f  >  <5  . 


(65d) 


Thus,  describes  that  unsteady  part  of  fluid  motion  which 
would  take  place  with  the  given  initial  conditions  and  a  body 
fixed  relative  to  the  translating  reference  frame. 

Lastly,  let  us  state  in  the  following  an  important 
property  of  the  integral  operator  defined  in  (57):  Let  fiC'f) 
and  'lf'(X)t)  be  any  integrable  functions  and  be  the  in¬ 
tegral  operator  defined  in  (57);  then 

To  prove  this,  let  the  following  substitutions  be  made 


in  (57): 
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(f,  (*;+)  =  (*;{-«)<!«  , 


etc . 

Then  from  (57)  we  have 

c[  ^«]  ^  4ir J ;  ^-wJ ^ 

+  II Js l  Ij'flii-H)  Gv(&’,t;o) 

So 

^(fjT  llJs  jT^^(l,r^(-!T-u) 

S0 

~]lJj  II  ^  t-u)^  [zfe  - 

~cj/r  ^'o,x>  T’i){£6(*:'?‘e’'*;t-T)j 

F 

~f  j/TtfP'^  >'r-u)  l 2  2*6+ 3  (2G)j 

+  J  J^Tj>  ^  J  T-«)j])G(&:tW,*‘T^ 

0  Pa  • 
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jr  l Tju  p(u  ]  /M  f~  t)  j s 

F  ~4^f(l>6)  J  Jf-Jj 

~Ci  Jr[  fr'0'* !T~H)  1^6(2 !  <’0,  r,  V-TjJ 

-  y y  <tx[  Juf}(u)(j)  ij'ff.o.  X  ;  T-n)j  z  ])ZG  +3C  yf  (bG)  j 

» y*  '*  Jr.  1 


Note  that  in  obtaining  the  last  expression  the  integral 
Jdu/}(n)  has  been  moved  outside  of  the  surface-  and  contour- 
integral  signs.  This  operation  is  justified  since  the  domains 
of  integration  ,  f  and  ft  are  independent  of  the  parameters 
T  and  W  .  Recall  the  following  formula  from  analysis,  some¬ 
times  called  Oirichlet's  formula:  If  /(r.  wj  is  continuous 
in  the  domain  O'-  T  £  :  ,  0 u  ■-  t  ,  then 
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where  w*  have  put  T  *  T-K. 

With  th«  use  of  this  formula,  the  last  expression  becomes 

1  ^ 

j'Ju/SOfi  4ir  if(£,  4-»)  Jh/H i-<*)  £v (a ;  2  ; o)  <h 

Sc 

+  J7-Jj  Lf(l  ;  T'J  6 x') 

■f  Jr'jj  ‘J'(-or,-r,)£lv1G(£:Uli:i'u'r')' 

-cj'Jup(*)J'  Jr’jj  ,,_f)  t'){2-C1(x^,o,x, 

F 

'  f  l  I  0  f  r)  \zj>Y,  +JC  £ 

•4  F, 

T  l^uftlu>J  /  if,  ~r')\])G 

Fc 

=  jf  <f  !  v} *•  u> , 


which  completes  the  proof  of  the  proposition. 
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V.  Decomposition  of  the  Velocity  Potential  Jf/  . 

Recall  that  in  (62)  we  wrote  down  an  integral  equation  to 
be  satisfied  by  which  can  also  be  written  as 


-  i  (  ^a(t)  jl  ^(t)O(s;  I;  V-t)  Js 

t>  ,  i 

*  t,k(±)  (?U;l:*-T)<k  .  (67a) 

s. 


Let  us  next  consider  in  the  following  two  boundary- value 
problems.  Suppose  that  (fFl  and  (f^'  are  two  potential 
functions  which,  in  addition  to  conditions  similar  to  those 
satisfied  by  the  potential  function  on  the  undisturbed 

free-surface  F,  on  the  bottom  S  and  at  infinity  (i.e.,  those 
stated  in  (35),  (36),  and  (38)  or  (38a)),  satisfy  the  following 
initial  and  boundary  conditions: 

I4n'(t;o)--‘0,  Iff  .  l’t-1)  <6,b> 


qF“’  (hi)  -  1  «*■' *->  ** "  •  <‘7c> 

Se  A  _  1 

fc/n(£i+)l  ®  for  +>°  ■  (67d) 

IS,  i*’ 
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By  precisely  the  same  analyst*  used  in  deriving  (58),  the 
integral  equation  for  the  function  ,  we  may  obtain  the 
following  two  equations  satisfied  by  and  (fn  , 

respectively: 

t‘>  °  c 

and 

l  tt<VG(i;l:i-T)Js.  <870 

Since  the  operator  q(  is  linear,  adding  (67e)  and  (67f)  together 
gives  the  equation 

i  j  tyf*1  +  tyF'*'\  (SH)  =  2  I  J t  Xi((ljjlit(i)6(*:m-r)‘ls 

k-'  *0  * 

^  (67g) 

jT6k(r)jJti(i>Cf(**:4  T)Js- 

So 

If,  as  we  have  assumed,  integral  equations  of  the  form 

have  unique  solutions,  then  we  can  conclude  that 


(67h) 
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.  u i  ft) 

Define  the  two  functions  <K^  and  <K^  to  be  the 
solutions  of  the  following  integral  equations,  respectively: 


and 


iW}(m) 


(68a) 


With  the  assumption  of  the  uniqueness  of  the  solutions  of  the 

Y/"’ 

equations  (6Ea)  and  (68b)  it  can  be  shown  easily  that  and 

'  satisfy,  respectively,  the  following  conditions: 


and 


(a 

K*  (*'■+> 


s. 


*  . 


for  ^  >e  j  (69a) 


k(!>  ■ 

4  *  /.  2. 


for  +  >  0 , 
4. 


(69b) 


Besides,  they  satisfy  also  the  same  free-surface  condition, 
conditions  on  the  bottom  and  at  infinity,  as  those  of 

We  now  claim  that  the  following  decompositions  of  the 

/  oi  /.  « > 

potential  functions  (ft,  and  lfFI  satisfy,  respectively,  the 
Integral  equations  (67e)  and  ( 67 f ) : 
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/  (II  &  /  t  ..  m 

%(*it)*LJ  °<i  h)  (i>  /-t)  Jr  , 


(70a) 


and 


i  /  f 


fy,  (*'<+)  I  (Z'I'T)Jt  ,  (70b) 

i‘i 

To  verify  this  statement,  let  us  apply  the  operator  o(  to  both 
sides  of,  say,  (70a).  Tl\en  we  have 


dffoi  =  2L  c(  j  l  -X((  t)X/(z;i-t)  <(s  j 

*  ^ l  *t(l]  Jr 

=  ^J^kh)  Jj  1 1^(1  )(,(&■  l; -I-t)  J$  , 

^ 0 

where  we  have  made  use  of  the  property  of  <£  stated  in  (66)  and 
equation  (66a).  From  (67h),  (7^8),  and  (70b)  it  is  clear  now 
that  the  function  has  the  following  decomposition: 

Ft  (-'•+)  -  i  /  J' 

k-i  * 

■+ 1  £  +-t)  JT ,  J 


(70c) 


We  now  show  that  the  decomposition  (70c)  toegether  with 
the  boundary  conditions  (69a)  ,  (69b)  and  those  satisfied  by 
<TU  and  d\ jj  on  the  undisturbed  free  surface,  the  bottom, 
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and  at  infinity  in  fact  is  essentially  that  of  Cummins  (1962) 
as  it  was  later  modified  by  Ogilvie  (1964).  Henceforth, 
reference  will  be  made  to  Ogilvie1 s  (1964)  paper. 

After  integrating  by  parts,  the  decomposition  (70c)  takes 
the  form 


On  the  other  hand,  Ogilvie  introduced  two  sets  of  functions 
and  proposed  the  following  decomposition: 

%(*■■*)  ■  jt  +  i  «,<*>%<» 

(ft,  i  [  (72) 

Thus,  the  decomposition  (71)  differs  from  that  of  Ogilvie 
in  only  inessential  ways.  Instead  of  starting  with  initial 
data  at  "t  =0  ,  Ogilvie' s  decomposition  starts  from  a  state  of 
rest  at  i* -oo.  In  fact,  the  functions  V/*  ^  »  VzV~^» 
Xifat)  ,  and  %**(*''*)  can  be  identified,  in  that 
order,  with  (*•*0)  >  anc* 

(&;i)  ,  respectively.  If  this  is  done,  then  according 

to  the  conditions  imposed  upon  the  functions  and 
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by  Ogilvie,  our  present  functions 


oust  satisfy  the  following  conditions: 


K01  and  K' 


o.  *?(' 

!  0  on 

f  ; 

(73a) 

XJ. 

n,tU) , 

Xt>~ 

f°lr 

; 

(73b) 

(h- 

*DC 

}K%  - 

0,  /»/. 

'2  ; 

(73c) 

K, 

>•  ; 

=  o  ,  (*;+) 

*  o  ; 

s. 

(73d) 

(&;  +  o  j  = 

i;  +o)  =  o  ; 

(73e) 

„ 

J'iit 

(x.o.1,,,+0) 

(n 

'0.  >,+<?)  = 

. 

(73f) 

Let  us  show  in  the  following  that  these  conditions  are 
indeed  satisfied.  We  see  that  (73b)  and  (73c)  are  already 
satisfied,  for  they  are  the  same  as  those  conditions  originally 
satisfied  by  <K^<"  and  in  (69a)  through  (69c).  The 

conditions  (73d)  follow  immediately  from  (69a)  and  (69b)  since 
both  fli  and  are  independent  of  time  here.  In  order  to 

verify  (73a),  we  need  the  equations  to  be  satisfied  by  <K*  (-  ~t0) 
and  (*  '<  *°)  ■  They  may  be  obtained  from  (68a)  and  (68b) 

by  setting  {  -  fO  .  Thus  from  (68a),  (68b),  and  from  the 


-53- 


operator  of  defined  in  (57)  we  have 

4r X* (it ; +o)  * H^dno) I  ;o)ds  ~  Ij ' tuk(i)G(is;l>oUs . 

£.  S0  (74a) 

4irXtl,(£;*o)  +  Ijxjl (1;*o)Gt,(*''1'°)Js  *  JJi^DGUH :o)Js. 

sa  s,  <74,b> 

Recall  that  In  (42b)  and  (43)  we  had  ;*)  *  O  and 

G(x;2;t)*s  so  that,  for  X  s  (n.o.^)  we  have 

6  7. r;  oj  =  GO-y.rtX’W’O)-  o  (75) 

for  all  1,  V,  J  .  But  then  also 

Gf(x°-y,  f. o)  -  y  for  -  on  f1  -(76^ 

Hence 

Gi-v^  0’^)  *?•  ?;  o)  =  G?  ■+  11  oi  C»;  -t  Ho?  -  o .  (77) 

Therefore,  for  *  =  on  f  ,  the  equations  for 

and  Xt‘(X;+0)  ,  i-e-,  <7*a)  and  (74b)  reduce  to 

and  ^(K.o.y,^o)^o  •  Thus  (73a)  is  also 
satisfied.  There  remain  the  last  two  conditions  to  be  verified, 
namely,  (73e)  and  (73f).  For  this  purpose,  we  need  equations 
to  be  satisfied  by  ^7*  - *0)  and  •  They  may  be 

obtained  by  differentiating  both  sides  of  the  equations  (68a) 
and  (68b)  with  respect  to  i  and  the  identity  (59).  For 
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instance,  the  aquation  for  ia  obtained  aa  followa: 

i. 

*  " 

F  1 
p f 

-  y  ^  ^  (l.o.s,+0)  \l)6(z;  t,o,\}f)  J  J ^  m 

But  we  have  Just  shown  that  to)*  O  (see  also  (73a)). 

Hence  the  last  equation  simplifies  to 

Or.*)  -  ///».,( l)6Vfi;£;«</i 

-  fj (it to)  Gvt (&>  i  /  f)  Js 

So 

*  c  II ^in(^°'  x>*°)  jif 

F  1 

Similarly,  for  we  have 


-55- 


(«;*)  -  ■f.tUs, 

S. 

„  ft)  .  >  (79) 

-  Hkk 

s, 

-+c Uni'll'0' f'+0)  ; 

f  ' 

For  /  s»  +0  ,  since  {r?  ^2£  ’  t/O.  T  ;  0j  -  G>t  (&  •  7'  ?  ' 

(78)  and  (79)  reduce  to 

(«;«)  +  \lk^(\:*o)Qv(^V.o)h^o  ,  (80) 

S0 

and 

H =  o.  (8i> 

So 

Obviously,  <*£  U  *0)  -  0  and  <k£  are 

solutions  of  (80)  and  (81),  respectively.  That  they  are  also 
the  only  solutions  is  shown  in  Appendix  II  of  Wehausen  (1965). 

The  last  condition  of  (73),  i.e.,  (73f),  follows  easily 
from  the  free- surface  boundary  condition.  Since  the  free- 
surface  condition  holds  for  all  i  ,  at  +  0  ,  in  particular, 

we  have 

Xttt  (*' *  h  K*lx  +  =  0 , 


56 


whara  1-  hi  and  ^  .  But  for  i • +o  t  wa  hava 

»  (*'°  y+0)=  0  «nd  (x>aj,,+o)-o  which  follow 

ianadiataly  from  (73a)  and  (73a),  raapactivaly.  Hanca,  from 
tha  fraa-aurfaca  condition  wa  hava 

°\kU  ( *'°\r.+o)  »  -  ,  A- 1.2  ;  i-  1.2.  . 

Thua  dCj1'  and  aatiafy  all  the  conditions  (73a) 

through  (73f). 
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VI.  The  Forcei  and  Moment!  Acting  on  the  Ship  Hull. 

In  the  following  we  shell  first  derive  expressions  for 
the  forces  end  moments  ecting  on  the  ship  hull  in  terms  of  the 
steedy  end  unsteady  velocity  potential,  and  (ft  .  The 

various  types  of  decomposition  which  were  made  previously  for 
the  potential  ty,  will  then  be  substituted  for  (f>t  . 

In  principle,  the  only  correct  way  of  calculating  the 
forces  and  moments  acting  on  the  ship  hull  is  to  integrate  the 
pressure  around  the  actual  wetted  surface.  However,  by  an 
approach  similar  to  that  used  in  obtaining  the  linearized 
boundary  condition  on  the  hull-fluid  interface,  it  is  found 
that  one  may,  in  practice,  integrate  the  pressure  around  the 
steady  reference  surface  S0  provided  a  proper  correction  is 
made. 

The  present  calculation  of  forces  and  moments  is  similar 
to  that  of  Appendix  A  of  Ogilvie  (1964)  in  the  following 
senses:  a)  the  same  assumption  is  made  about  the  wall- 
sidedness  of  the  ship  hull  in  the  vicinity  of  the  equilibrium 
water  line  where  the  undisturbed  free- surface  intersects  with 
the  undisturbed  ship  hull;  b)  the  same  criterion  is  used  in 
discarding  terms  of  higher  order.  An  effort  will  also  be  made 
to  preserve  the  same  use  of  notation  so  that  any  cross 
reference  between  the  present  work  and  that  of  Ogilvie  (1964) 
may  be  easily  made. 

To  begin  with,  we  shall  consider  simultaneously  the  body 
reference  frame  and  the  steadily  translating 
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reference  frame  Oxjy  .  Toward  the  end,  however,  only  the 
Oxjfj.  -  system  will  be  needed. 

Ae  before,  let  S:  F(*>1$)  =  0  ,  in  the  body  reference 
frame  6t$$,  ,  describe  the  actual  ship  surface  which  is  below 
the  equilibrium  water  line.  Then  the  equation  F  (*'}•}■)*  & 
in  the  steadily  translating  reference  frame  Oxyjy  describes 
the  reference  surface  which  coincides  exactly  with  S 
when  the  latter  is  in  its  undisturbed  position.  Note  that 
both  S  and  S„  are  defined  by  the  same  function  F  but  in 
two  different  reference  frames. 

For  convenience,  we  shall  define  the  following  term, 
a  point  Qo  is  the  'image  point'  of  a  point  Q  on  the  ship 
surface  S  if  the  coordinates  of  in  the  0  -frame 

have  precisely  the  same  values  as  those  of  Q  in  the  0x}£ 
frame.  Thus  Q  coincides  precisely  with  Q„  when  S  coincides 
with  S„  •  Suppose  that  Q  is  the  unit  inward  normal  to  the 
ship  surface  S  and 

fj.  ( =  n*  is  ~  F>;  > 

A,  i  =  I.1 ,  3. 

If  H o  is  the  unit  inward  normal  to  S„  ,  then  due  to  the 
way  S»  is  related  to  S  we  have 


n.;  (*•},})  5  /?.-«-  =  > 

A,  *  »  A  2,  3. 
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Thu*  hi  and  tin  are  in  fact  given  by  the  aame  function. 
Hence  if  QiS  has  the  coordinates  (a.b.c)  in  the  0 
frame  and  Qo  £  S,  is  the  image  point  of  Q  ,  i-e.,  Q  has 
the  same  coordinates  (a.b.C)  but  in  a  different  reference 
frame,  Oxjy  ,  then,  clearly, 

f)i(Q)  -  ,  (82) 

i  =  /•  2.3. 

The  force  on  the  ship  hull  is  given  by 

X  =  j] I’ll  Js  ,  <“> 

A 

and  the  moment  with  respect  to  the  point  0  ,  the  origin  of 

the  body  reference  frame  is  given  by 

1JL  =  //  f’X'xn  Js  <e4) 

where  X"  =  X  &  *  $  £i  ■+  £  £  is  the  position  vector  of  a  point 
on  the  BhiP  hul1  with  reference  to  the  body 
reference  frame  Ox$£  and  Sw  is  the  actual  wetted  surface 
of  the  ship  hull. 

In  order  to  use  the  results  of  the  present  calculation  of 
force  and  moment  to  write  down  equations  of  motion,  we  shall 
resolve  the  force  and  moment  along  the  steadily  translating 
reference  frame  0x^  .  For  convenience,  let  us  consider 
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fir»t  the  decomposition  of  the  force  vector  along.  From  (43) 

we  have 


X>ej  *  jj  f>  h 


3  tv 


- lit  lift -  Ji 


Sy 


if  i $ 

s 

+  tit  I 

s, 


(85) 


where  *  S  U Si  ,  and  S,  *  ^  S. 

Note  that  in  obtaining  (85)  we  have  used  the  result  from  (5), 

namely,  §i  •  §j  *  Stj  •+  . 

Let  Q,l  S,  be  the  image  point  of  Q  l  S  -  Suppose  that 

Q0i  Ss  has  the  coordinates  (*“}'})  in  the  @X}J> -frame; 
then  <3  has  the  same  coordinates  (*'}'&)  in  the  -frune. 

The  coordinates  of  Q  l  S  in  the  0* U  -frame  are  then  given 
according  to  the  transformation  (3)  as  Q  =  (  *.•  +  */*•+ X*  )  . 
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We  8 hall  assume  that  the  function  can  be  expanded  into 

Taylor's  series  as  follows: 

*  t  (*■  M  +  l*'i 4  ^ . *%*»]  (*'l'3"f)  +  CM)  , 

which  means  precisely  that 

/f>  (Q;U=  /  fait)  -t  [Xu  +tXtnn  +  (86) 

Thus  from  (86)  and  (82)  we  may  write  the  integrand  of  (85) 
as  follows: 

V  .  >  (87) 

=  j  hoj(Q.)  +  Ljt;  $h  lh;(M  1  ;  +)  4 

■»  Hj  (Q,)[*u  +  £/*.„  $rn  Xn  ]  #(£»)- 

Although  the  expression  obtained  on  the  right-hand  side  of  the 
equality  in  (87)  is  a  more  complicated  one,  it  has  the  advan¬ 
tage  of  being  evaluated  on  a  prescribed  reference  surface  S. 
which  is  stationary  with  respect  to  the  Ox^-frame.  It 
follows  immediately  from  (87)  that 
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IJ I  n.  (Q)  t  fy;  A*  k  (Q)}  f  (Q-  i)  h 

s  * 

-  jj  {  n,j(4.)  +  ^  \  ki  J  j (*•:  t)  <k 

S'. 

* // *u  4  t-lmn  Ah  x-ltjfaiyJs  +OUii), 

o 

which  shows  how  an  integral  over  the  surface  5  may  be 
approximated  by  Integrals  over  the  stationary  (relative  to 
the  frame)  reference  surface  S,  • 

We  need  next  an  appropriate  expression  for  /  •  From 
Bernoulli1 8  equation  we  have 

{{*§)*•  const., 

where  we  let 

-ex*  ft  OD-t  (?,(*•+).  (89) 

Assuming  h  0  .  7#  “  C  and  0  at  *  =  4  00  . 

|s  o  ,  we  have 
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f(W'i.t)  =  -fn~  i  t( v$)z  4  i 

We  shall  follow  the  usual  practice  of  linearized  theory  of 
discarding  the  last  two  quadratic  terms  in  (90).  However, 
one  should  realize  that  omission  of  the  term  -£■  f  (V(fof  may 
not  be  proper  for  the  case  of  a  deeply  submerged  body.  Let 
us  then  write 

t(i-’f)  =  -m  ■>  fc(f,x  - 

(91) 

This  is  now  to  be  substituted  into  (88)  in  place  of  /*  .  If 
this  is  done,  then  we  have  from  (88)  and  (91)  that 

s  * 

-  fjj  j-n[noj 

s, 

-rivKjrcfaWr.t)*  na)-H] 

+  (U  4  S (x,*h . 
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In  (92)  we  have  discarded  terms  like  fi  l  ftt 

*  F^/], etc- »  which  are  of  the  order  0(£m)- 
Next,  let  us  consider  the  integral  over  S,  in  (85). 
Similar  to  the  relationship  between  $,  and  S  ,  let  So,  be 
the  area  on  the  undisturbed  position  of  the  ship  hull  which 
has  the  same  size  and  form  as  Sf  .  The  area  Se,  can  be 
determined  from  S,  in  the  same  manner  as  finding  the  point 
Q,  on  S,  after  a  point  Q  on  S  is  assigned.  The  free- 
surface  elevation  }  -  YL*'b*)  can  be  expressed  in  terms  of 
the  ship  coordinate  system  0x‘^£  as  follows: 


■+  h.o.t., 


where  we  have  used  the  transformation  (3)  and  expanded  the 
function  X  into  Taylor's  series.  Thus  from  the  last  expression 


we  have 


Y(t'lt)- +  +  j]Y*  4 

r  ,  .  i  >  (93a! 

+  +  Yj,  +  h.o.t. 

sy 

Thus,  in  the  OSjtf  system  S,  is  defined  as  that  part  of  the 

hull  surface  ?(*)•?)-- °  bounded  by  the  equilibrium  water 
line  j  =  0  and  the  free-surface  elevation  %=  YCx'l’^) 
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according  to  (93a).  Hence,  in  the  O*"}}- -system,  Set  is 

‘i* 

obtained  as  that  part  of  the  reference  surface  F  *  0 

bounded  by  the  undisturbed  free-surface  }=  0  and  the  wavy 
surface 

;  *  9  (*•}■  +  ) 

(93b) 

The  relationship  (87)  permits  us ,  as  in  (88) ,  to  replace  the 
integral  over  S,  by  integrals  over  S„  as  follows: 

II  f  (M)  /  ^  (V  +  fy*)}  <*S 

-  /j  t  (Qo.i)  I  Js 

S0( 

+ IJ *  oc^) 

S‘l 

For  simplicity,  let  us  further  assume  that  the  ship  hull 
is  wall-sided  near  the  equilibrium  water  line  so  that  the 
integrals  over  Sot  in  (94)  can  be  written  as  the  following 
iterated  integrals.  The  first  integral  on  the  right-hand  side 
of  the  equality  in  (94)  then  becomes 
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Uft*W)  I  HfC**})*  LikiW.i  |  ^ 

-  J  J/i\ fa }(  f  (*’ V>i"t)  J} 

y 

-/  Jj{n.i(*.°+)*£^6tii..:}J'  f-tn  *  f‘b,(*w- 

where  To  is  the  contour  where  the  undisturbed  free-surface 

A 

intersects  with  the  reference  surface  S„  >  and  Y  is  8iven 
by  (93b).  Also,  in  (95)  we  have  used  the  expression  for  / 
from  (91)  and  set  (*'}’£)  s  due  to  the  wall- 

sidedness  of  the  ship  hull  near  the  equilibrium  waterline. 

The  integration  with  respect  to  the  variable  Jf  in  the  last 
expression  of  (95)  may  be  carried  out  as  follows: 

A 

l  hn+fc L (K id -f[(7rc&)9,( *w) + vi -vfi]j 
=  -in  fz(*'ht)  +  fc  ?  ij* '/<?']')  - 

-f  ?■[&-<&) 9,  (*-/*'?■>. *)]■*  0(0  >  <96a> 

where  0  yt/ £  l 

Since 
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we  have  from  (93b)  that 

(96b) 

+  f  [  *, +  A  J-  W  7  ( *'»l  )  +  f  [*•+*•}- 1*1  by  *  0(E" }- 

We  may  further  expand  those  functions,  (foK (*’/* Y>  , 
tyr (X'/V-  »  and  fe  <V?' etc->  lnt0  Taylor’s  series 

as  follows: 

* 

if,  (*./?.  i.t)  -  «w*.**f) 0(ti)> 

etc. 

Substituting  (96c)  and  (96b)  into  (56a)  we  have 

=  'f  fC  %,(*•*' })[(£'*&)  ^  *  • 

This  result  is  now  to  be  substituted  back  into  the  last 
expression  of  (95).  Then  we  have 

II H*'  lb *)  1  n.j  (*>H)  +  ty; J 5 

Set 

=  -  fc  I  +0(0. 

To  ?  ^  (97) 
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Following  the  same  procedure  as  above,  one  may  obtain  the 
following  reault  from  the  eecond  integral  in  (94): 


//  tl.j  (*')•})  LtmnLKlt.l  (*>>'}■+)  JS 
S.< 

-  -fcj  Wt*-Ab]  U  +  OUZ) .  (98) 

The  appropriate  firet-order  expression  for  the  force 
component  is  now  obtained  by  adding  together  (98),  (97),  and 


(92)  as  follows: 

1:  =  f  JJf-Ulty  +  fyA  *'**'*' *'*]' 

*  S. 

-  ruf  [  & - « 4)  9,  (*'Ui)  +  H(*'U)  ^9,] 

+  C[('Uj+tyify'1«)hx  f  JS 


V  (99) 


-  tef  n0j (m)  (t?  ~c&)  #  (*'0>h*)+ 


Note  that  (99)  corresponds  to  the  expression  for  on  p.  99 
of  Ogilvie  (1964)  except  that  -Xj  in  (99)  is  now  resolved 
directly  along  the  steadily  translating  reference  frame  0*n  . 
For  interpretation  of  terms  appearing  in  (99)  and  remarks  about 
the  correction  term,  i.e.,  the  contour  integral  around  T,  in 
(99),  we  refer  to  Ogilvie’s  (1964)  discussion  following  imme¬ 
diately  his  expression  for  •  However,  we  call  attention 
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here  to  the  fact  that  the  correction  term  is  of  order  not  lower 
than  £s  &n  and  vanishes  entirely  when  the  body,  is  completely 
submerged . 

We  may  proceed  in  the  same  manner  to  obtain  an  expression 
for  the  moment  components  %(;  .  As  before,  let  Q  be  a  typical 

point  on  S  and  Q0  be  its  image  point  on  .  Define 

L  =  x*n„  =  n.j+3  £■  ,  (ioo) 

—  /«/  f  * 

where  %  -  *  ^  1  4  £  -■*  is  the  position  vector  of  the 

point  Qo  £  So  in  the  0*^- frame.  Then  the  position  vector 
of  the  image  point  C(  i  S  in  the  Oxff- frame  is  given  by 
X'  =  V  £)  +  y  }  £s  .  Let  us  define  also  that 

A  =  x'  xn  =  t  h+3  ii  .  (ioi) 

* 

Thus  from  (101),  (100)  and  the  fact  that  Qo  l  is  the  image 
point  of  Q  2  S  we  may  again  establish  the  relationship 

n^(Q)  =  nci+3iQo) ,  /  =  i.2.3.  do2) 

We  shall  again  decompose  the  moment  acting  on  the  ship  hull 
about  the  origin  of  the  ship  coordinate  system  Oxfp  directly 
along  the  steadily  translating  frame  as  follows: 
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7rij  *  7n- S{  »  jj  f  s}-  <4 

-  Js 

-  H 4 1 4  £><-  s»  ) Js 

Sty 


*•}>  k  -  I ■  *>  •*• 

Note  that  (103)  is  essentially  the  same  expression  as  (85) 
except  that  the  index  j  is  now  replaced  by  J+  3.  With  the 
relationship  established  in  (102),  namely,  Dpi  (Q )  ~ 
where  Q  f  S  and  Q,  f  S.  ,  the  same  analysis  used  in  obtaining 
the  final  expression  for  Xj  holds  also  for  .  For 

convenience  let  us  write  Mj  * -fy#  '•  then  <103)  wil1  lead  to 
an  expression  precisely  the  same  as  (99)  except  that  Xj  , 

,  and  //,,  are  now  replaced  by  ,  ^»i*3  *  and 

t l,^i  .  Thus  from  (103)  we  have 


=  fll  h  f  dvr’  +£>  *  1  L:«  1 + ^ *”  ^ 

s» 


(104) 
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+ c  [(rkp-ttyiktuHi)  fc  +  )ktPs 

-fc j  n.j+3  (*> °-P)  fo* ( *■ °4)  \jQi~Cih)(f,(x'0-l>t)+ 

*  y>  +  *~  &  i  ]  J*, 

*,  j>  k  *  A  2,  J . 

Note  that  from  (100)  we  have 

h,^(*-  ll)  -  fyji  **  n.;(*W)  ,  J.j.  *  *  1,2.3.  (105) 

Since  f1*j}(X'].l)athj(*‘0'],),.‘*i.2.3,  follows  from  the  wall’-sldedness 
of  the  ship  hull, and  in  the  region  S,,  we  have  111  -  I)' (*■}■+)!  , 
the  difference  between  (*'}•))  and  /7.y+j  (*•»•})  is  of  order 

<Ss  or  £M  .  Thus  the  contour  integral  in  (104)  is  the  only 
contribution  to  the  first-order  expression  for  the  moment  com¬ 
ponent  from  the  correction  term,  i.e.,  the  integral  over  S»i 
which  is  similar  to  (94).  Furthermore,  the  wall-sldedness  of 
the  ship  hull  implies  that  lr)m  o  for  £•/ 1  so 

that  °'l)  -  Hot  (*'  °'l)  -  0.  Thus  the  contour  integral 

in  (104)  vanishes  except  for  ^  -  2.  Hence,  it  contributes  to 
the  yawing  moment  alone. 
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Lastly  let  us  substitute  into  the  expressions  (92)  and 
(104)  the  various  types  of  decomposition  we  made  for  , 
namely. 


ban) 


*  h(l;i)+ljm  ZtMXl'U '•*’*)  Jr 

+  '(*;  t-r)  Jr  .  ! 


(106) 


Note  that  to  be  consistent  with  the  use  of  the  notation  cX  for 
the  displacements,  we  should  also  write  Q+  *  °(^3  »  1>2»3» 

and  Jl,  a  o(2  • 

With  the  substitution  for  ^  J-rom  and  use  of  the 

boundary  condition  fl,j  (l)  ~  Xtn(£.:*o)  for  £  i  S„  from 
(69a),  the  equations  (99)  and  (104)  lead  to  the  same  equations 
in  the  following  form: 


J.  =  lj0  -i  | 

-Z  (t)  Ljt  U- -r)  Jt  -  (r)  Mjt  ((- tW"  , 


(107) 


where  i  ■  1,  2,...,  6, 


s. 


(108a] 


(108b; 
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s. 


MjtH)  -  f •Jj^(^-o)\(h-c»)^i^  *  v<f-‘  7i(t'\ Js 

Si 

*  j  fcj  ck)X+ 

P 0 

* 

Cft  =  fl/k^  (**• °>lu  «•  -fiiwj-  c(®>  •t#.,}  * 

+  fcJ^ k^(*’0,},4°)\[e>’§*(u] tyox (Xl 0t ^ 


where  <a^  [  £«  *  §i  (*)]  ~  $>  + 

i- 1 

§t ,  i-  1-2, 3, 

£i  -  I 

l  ^/i-3  *  >  $  m  4i  S 1  £  , 


(108c) 


V  (108d) 


l  C108e) 


(108f) 
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It  may  be  worth  mentioning  again  that  -2.^'  given  by  (107) 
represents  the  force  and  moment  (about  O ,  the  origin  of  the  ship 
coordinate  system)  acting  on  the  instantaneous  position  of  the 
ship.  One  can  see  from  the  definitions  for  the  various  quan¬ 
tities  given  in  (108a)  through  (1081)  that  X.^0  represents 
steady  force  and  moment  components,  !■*  a  constant 

depending  only  on  ship  geometry,  and  £/*  are  constants 
depending  upon  the  ship  geometry  and  the  constant  average 
forward  speed,  and  Lji  and  Mjk  are  functions  of  time,  geometry 
and  the  constant  average  forward  speed.  None  of  these  quantities 
depends  upon  the  unsteady  oscillatory  motion  of  the  ship.  (We 
note  that  bjk  ,  LJk  ,  Htt  .  and  Cjk  do  not  have  the  same 
meanings  here  as  in  Ogilvie's  paper.)  includes:  a)  the 

static  buoyancy  and  wave  resistance  of  the  ship  moving  with  a 
constant  forward  speed  in  its  equilibrium  position;  b)  the  total 
hydrostatic  and  hydrodynamic  force  and  moment  acting  on  the 
ship  due  to  its  own  motion;  and  c)  the  force  and  moment  resulting 
from  the  action  of  the  incident  and  diffracted  waves  upon  the 
ship. 

-  It 

We  remark  in  passing  that  the  quantities  Ojt  defined  by 
(108g)  arise  only  if  the  force  and  moment  vectors  are  resolved 
along  the  steadily  translating  reference  frame  Ox }}  .  Thus, 

if  we  set  CjfO  for  j-i  -  1,  2 . 6,  then  the  expression 

(107)  will  give  the  components  of  force  and  moment  resolved 
along  the  ship  coordinate  system  0  *  . 


-76- 


VII.  The  Equations  of  Motion. 


Without  the  loss  of  generality,  we  shall  suppose  that  the 
center  of  gravity  of  the  ship  is  located  at  a  point  with  the 
coordinates  (  0  ,  %  ,  0  )  referred  to  the  ship  coordinate 
system  0  >'■/£.  .  We  shall  also  assume  that  the  propeller 

thrust  T  i®  directed  parallel  to  Ox  along  a  line  of 
action  which  is  tr  units  below  the  center  of  gravity  and 
lies  in  the  -plane. 

In  writing  down  the  equations  of  motion  for  the  ship  the 
following  result  of  Euler  in  rigid-body  dynamics  may  be  used: 

The  motion  of  a  rigid  body  can  be  determined  by  treating 
separately  the  motion  of  the  center  of  gravity  and  the  rotations 
motion  about  the  center  of  gravity.  Moreover,  the  motion  of 
the  center  of  gravity  can  be  determined  as  the  motion  of  a 
particle  of  the  total  mass  // 7  of  the  body  subjected  to  the 
total  force  [_  applied  to  the  body,  and  the  rotational  motion 
of  the  body  about  its  center  of  gravity  can  be  determined  as 
if  this  point  were  fixed  and  the  body  were  subject  to  the 
moment  of  the  applied  force  about  this  point.  If  we  use  the 
steadily  translating  coordinate  system  Oxyj  as  our  inertial 
reference  frame,  then  the  motion  of  the  center  of  gravity  is 
given  by 


m  *,  =  F, 
m  j t,  —  Fi 


(109! 


where  £  -  Fi  €1  +  Ft  §*  i  F»  §1  is  the  totel  external  force 

■  I  M 

acting  on  the  ship  and  X,  ,  J,  ,  and  ^  are  the  three  com¬ 
ponents  of  the  rectilinear  acceleration  of  the  center  of 
gravity  in  the  0*  0y-t  and  -directions,  respectively.  In 
our  present  generalized  notation  we  may  write  these  quantities 
as  X,  *  <*,  ,  jf*  =  ,  and  £  =  o(3  .  The  motion  of  the 

ship  about  the  center  of  gravity  is,  then  determined  by  Euler's 
equations  of  motion: 


I,  ~  J-n  ~  1)3  £<Jj  *f  Wj  [~-?/3  bJ/  ~  In  *  %s  £<Jj  ]~ 

oN 

-  6Jj  [~  In  &J,  4  1,  (*Ji  ~  It }  k)}  j  *  7//$  i  > 

“I/tUJi  ~  Iz3^0j  4  Wj  f  1/  kJi  "  I-  >t  ^2  — 

aT' 

-  [~l,j  U),  -  J/6  COl  -i  Jj  i‘Ji  J  —  'W$z 


/3  bJi  ~  1 23  LO  l  "t  I 3  \  III  *  1 1  ~  1 13  ^jj~ 


where 


-LOi[l,u),  -  Ju  <0,  -  ln  c03  ]  =  lih 


OS 


l  (110) 


K  =  Ai  $  +As  l 


(in) 


is  the  moment  of  the  total  external  force  £  about  the  center 
of  gravity  and 
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is  the  angular  velocity  of  the  ship  as  a  rotating  rigid  body. 
In  order  to  define  the  quantities  ,  let  us  introduce 

the  function  »  c^e  density  of  the  material  of  the 

ship.  Then  evidently 


M  =  UJ  fs(*-W)Jv 

V 

and  the  Ij/s  are  given  by 


i  -III ^  ~lllt‘lv+~*'Vv  • 

V  u 


y  aw) 


4  •mw-w 

V 


hi  -llJf.ixJv. 


where  V  is  the  total  volume  of  the  ship,  including  the  part 
above  the  waterplane. 

Let  us  now  show  that  the  quantities  ^0*  may  be  identi¬ 
fied,  to  the  first  order  of  approximation,  with  the  quantities 
Gj.  which  we  have  introduced  in  the  previous  chapters.  It  is 
a  well-known  result  in  the  kinematics  of  a  rotating  reference 
frame  that 


JeJJt  =  Sjji  to*  if 


(114) 


for  an  arbitrary  rotational  motion  of  the  coordinate  system 
dSyJ.  .  On  the  other  hand,  for  small  angular  displacements 
of  the  coordinate  system  we  have  the  relationships 


£*  *  £-■  +  ■+  0 (&•**) 

and 

-  %  +  ty.hL  ->  o(z'm) 

which  can  be  derived  easily  from  the  transformation  (4). 

Thus  we  have 

Jgi/Jt  •  t.ji &  e,  -i-OM)  - 

•  tiAij  -o<c>,'  (iw 

where  we  have  assumed  that  the  $kS  are  also  of  the  order 
.  Hence  from  (114)  and  (115)  we  have 


bJ/c.  -  (9k  +  O(c-m)  ~  ^11+3  ^ 

i  =  J. 

Let  us  now  write  the  equations  of  motion  for  the  center 
of  gravity  of  the  ship  according  to  (109).  To  obtain  the 
total  external  force  acting  on  the  ship,  we  must  add  to  Xj 
the  other  forces  acting  on  the  ship,  for  example,  the  gravity 
force,  propeller  thrust,  forces  due  to  wind  gusts  and  arti¬ 
ficial  restraints,  etc.  Thus  the  equations  for  the  ship 
corresponding  to  (109)  may  be  written  as 


-80- 

wh*re  Tj  represents  the  three  components  of  the  propeller 
thrust  resolved  along  the  0*-t  °f~t  an<*  -directions, 
respectively,  and  by  our  assumption  about  the  line  of  action 
of  the  propeller  thrust  we  have 

1 5  =  I-C;  =  (118) 

/*  J’ 

C-(t)  represents  all  other  external  forces  besides  the  force 
Z-(i)  ,  the  propeller  thrust  T-  ,  and  the  gravity  force 
-mL.  The  expression  (117)  can  oe  further  put  into  the 
following  form  if  (107)  is  substituted  for  X}«)  : 

Z  ("tji  «*  (t)  +L  +  h  ci*  ** cn 

*g/  «4(r>LjtU-T)</ T  l  ^'Tj</T 

=  +  ~Tj  (+)  —  My  -  Xiv-  (t)  , 

i*  '■  z>  3. 

where 

/n  o  o  o  o  o\ 

0  ho  0  0  o  o 

0  0  0  O  o)  , 

;-/•  a.i  ;  4  =  '• -•  -•  6. 
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We  shall  next  write  Euler's  equations  for  the  rotational 
motion  of  the  ship  about  its  center  of  gravity  according  to 
(110).  Let  be  the  vector  directed  from  the  center  of 

gravity  of  the  ship  to  a  point  <5  on  t*1®  wetted  ship  hull. 
Then  the  moment  about  the  center  of  gravity  of  the  ship  due 
to  the  pressure  distribution  around  the  wetted  ship  hull  is 
given  by 

Ij  +  r*  xn  Js  . 


But  by  the  assumption  about  the  location  of  the  center  of 


gravity  we  have 


K  -  x  -  %  &  , 


(121) 


where,  as  before,  2'  is  the  position  vector  of  the  point 
on  a  ship  hull.  Hence 


jj  j  tfi  x  tl  </i  =  Jj  £'  >  JJ  Js  -  jj  t  £*  ^ 

Sm/ 

=  j  2j„  ^  I  Z„  }  £j  ,  (122) 

"»  *  /,  Z,jt 

where  the  repeated  indices  imply  summation  and 
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Clearly,  Xj+s  and  -£j.  are,  respectively,  the  components  of 
the  moment  (about  0  )  and  the  force  due  to  the  pressure  dis¬ 
tribution  when  those  two  vectors  are  resolved  along  the  ship 
coordinate  system.  As  it  was  remarked  at  the  end  of  the  last 
chapter  that  the  and  the  may  still  be  given 

by  the  expression  (107)  if  we  put  C , ^  s  0  ,  &  m  1,  2,.  ..,  6 

in  (108g) .  This  is,  in  fact,  equivalent  to  replacing 
by  Cjk  in  the  expression  (107).  Let  us  denote  by  TJ+jWthe 
three  components  of  the  moment  about  the  center  of  gravity  due 
to  the  propeller  thrust  T*  and  by  that  due  to  all  the 

other  external  forces  besides  the  pressure  distribution  and 
the  thrust.  Then,  by  neglecting  higher-order  terms,  we  may 
write  the  equations  corresponding  to  (110)  in  the  following 
form: 

Or,  with  the  use  of  (107)  this  last  equation  may  be  further 
written  in  a  form  parallel  to  (119): 
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and 

"7^+3  -  ( 0  •  0  <  T  /r)  ,  1  =  h  *'  i  • 

Note  that  (119)  and  (123)  together  form  a  set  of  six 
integro-differential  equations  for  #(+(+)  ,  km  1.  2,...,  6. 
Such  equations  are  generally  solved  by  method  of  the  Laplace 
transforms.  Thus  the  information  about  the  position  and 


84 


velocity  of  both  ihlp  end  fluid  it  the  initio  1  instant  t-  o  , 
i.e. ,  oik  (a)  ,  0(4(0)  »  Y (*>l<o)  .  *nd  Yt(*^°)  »  t08ether 

with  the  Integral  aquations  for  fal*;*),  X*  (*•  *)  »  and 
X*  (*'>+)  ,  end  the  above  set  of  lntegro-differential 

equations  uniquely  determine  the  behavior  of  both  ship  and 
fluid  at  later  instants  of  time,  provided  one  grants 
uniqueness  of  solution  of  the  integral  equations. 
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